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TECHNICAL MEMORANDUM 


ON LONG-TERM PERIODICITIES IN THE SUNSPOT RECORD 


I. INTRODUCTION 


In recent years, substantial interest has been given to the sunspot record, 
especially its historical aspects. Much of this work sought to associate possible 
climatic excursions with variations in sunspot number. For example, Eddy [ 1—8] has 
done considerable research on the historical record of sunspot data and solar activity 
(also [9-43, 162-167, 172-183]). His studies, based primarily on naked-eye and early 
telescopic sunspot observations, fossil radiocarbon levels in tree rings, and auroral 
reports, have suggested possibly strong climatic associations in the Sun-Earth record; 
in particular, he has found evidence for the occurrence of prolonged "sunspot 
minima," e.g., the Sporer Minimum (A.D. 1400-1510) and Maunder Minimum (A.D. 
1645-1715), and for the occurrence of a prolonged "sunspot maximum," i.e., the 
so-called Medieval Maximum or Grand Maximum (ca. , A.D. 1120-1280). In addition, 
he has given tutorial remarks concerning the sunspot record post-1850 and pre-1850; 
i.e., he notes that the highest or best-quality sunspot data are those dated from 
about 1850 to the present, since it was about 1850 when solar observers first began 
to record telescopic sunspot observations in a systematic fashion, having just become 
aware of the existence of the cyclic nature of sunspots (Wilson [44] and Wilson et al. 
[45]). Sunspot observations prior to about 1850 are viewed as being somewhat less 
reliable, or at least somewhat different, as a simple comparison of selected cycle 
parameters for pre-1850 to those post-1850 will show. (Sonett [36] has suggested 
that there may be "bookkeeping" type errors in the compilations associated with these 
early sunspot cycles; also it must not be forgotten that daily sunspot records were 
not maintained on a continuous basis until January 1849, Waldmeier [170].) Cycles 
have been arbitrarily assigned numbers back in time into the 1700's (and before) , 
with sunspot cycle number (SCN) 1 having a maximum in June 1761; the present 
cycle, now in decline, is denoted SCN 21. 

Another aspect of this work has been to estimate or "predict" sunspot numbers 
as a function of time for future epochs, based on modeling or "post-dieting" of past 
cycles. This aspect usually incorporates various time series techniques which some- 
times make use of power-spectrum analysis results. Power-spectrum analyses suggest 
the existence of several possible short-term and long-term periodicities in the solar 
cycles; e.g., several authors have reported the presence of 5-, 8-, 10- , 22-, 90- , 
and 180-year periods, as well as the well-known 11-year average period cycle dura- 
tion in the sunspot record. Some authors have hypothesized that these periods may 
be due to planetary motion "tidal" effects, while others argue that sunspot-like cycles 
can be simulated with random noise. It must be pointed out, however, that some 
skepticism remains regarding the authenticity of most of these short-term and long- 
term periodicities, and, following Wallenhorst [62], it is important to note that the 
11-year periodicity is the only statistically significant one present in the sunspot 
data. We should recall that sunspots have been routinely and systematically monitored 
for a mere 130 years or so; thus, few or no cycles of periods substantially longer 
than 11 years have been observed. Wilson [44] , Wallenhorst [46], Newkirk and 
Frazier [47], Bray and Loughhead [48], Noyes [49] , and Meadows [50] have given 
additional comments regarding sunspot numbers, their history, measurement, and use. 



A number of references [51-161, 170-183] are included which allows the reader to 
follow the past 80 years growth in our understanding of sunspot cyclicity and to see 
the methods employed for estimating contemporary values of sunspot number. 

The purpose of this report is threefold: (1) to compare observed sunspot 

numbers at cycle maximum of contemporary cycles (i.e. , cycles 8 through 20) with 
similar computed values, based on the assumption that sunspots are periodic in nature 
(i.e., sinusoidal); (2) to extrapolate the preferred, empirical periodic function "back- 
wards" in time to qualitatively compare computed values with values estimated by 
others for such epochs as the Maunder Minimum; and (3) to extrapolate the preferred 
periodic function "forwards" in time to estimate smoothed sunspot number at sunspot 
maximum (denoted, for cycles 21, 22, and beyond. In addition, using linear 

relationships identified in Wilson [44], a number of cycle-related parametric values 
will be estimated and compared with observed values (when they are known) for cycle 
21 and will be projected for cycle 22. In Section II, the approach used in this study 
is described; in Section III, the results and discussion of the findings are presented; 
and in Section IV, the conclusions are stated. This report is a continuation of 
studies investigating solar activity and its estimate (e.g. , Wilson [44, 102], Wilson 
et al. [45, 100], and Teuber et al. [107]). 


II. APPROACH 


Since Eddy [4] (and before him McNish and Lincoln [123]) has cautioned that 
cycles observed prior to about 1850, corresponding to SCN 9, may not be reliably 
known, many solar statisticians today limit their remarks concerning cycle-related 
parameters, such as jascent period (dubjbed ASC , the time in months measured from 
R MIN occuri>ence 1:0 r mAX occurrence ; r min * s the minimum smoothed sunspot number 

or value marking the beginning of a cycle and is the maximum R^ value 

in the same cycle), descent period (dubbed DES, the time in months from 

occurrence to subsequent cycle Rjyjj jq occurrence), cycle duration (denoted MIN-MIN 

PERIOD or the cycle minimum-to-minimum period, numerically equal to ASC + DES), 
and others to cycles 8 through 20. Table 1 lists values for the cycle- related 
parameters for cycles 8 through 20. It also gives mean and "uncorrected" standard 
deviation for each parameter. Definitions of these parameters are given in the list 
of selected abbreviations found at the beginning of this report. 

Using R max values for cycles 8 through 20, this report compares a linear- 

regression approach and modifications to it (APPROACH 1) with two cyclic approaches 
and modifications to them (APPROACH 2 and APPROACH 3) . The two cyclic 
approaches presume a 50-cycle or approximately 550-year periodicity and a 90-cycle 
or 990-year periodicity in sunspot number, respectively. The reasoning behind the 
selection of these two periods will become evident in the following. 

If restricted to cycles 8 through 20, it is observed that the mean value for 
Rmax * s H6.2 with a 1-sigma confidence value (or standard deviation, s) of 36.7. 

(For small sample size, N, when N < 50, the "corrected" 1-sigma value is equal to 

/N m - 1) x s, or 38.2.) Thus, any given cycle can "crudely" be predicted to 
have, assuming a 1-sigma confidence with no apparent trend in the data, an Rjyj^ x 
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TABLE 1. SELECTED VALUES FOR CYCLES 8 THROUGH 20 


MIN— MIN 


SCN 

r MAX 

r min 

r mean 

S MEAN 

r chm 

asc chm 

d chm 

ASC 

DES 

PERIOD 

slope asc 

SLOPEdes 

8 

146.9 

7.3 

66.9 

43.4 

77.1 

23 

53 

40 

76 

116 

3.490 

-1.795 

9 

132.0 

10.5 

55.9 

30.4 

71.3 

44 

39 

55 

94 

149 

2.209 

-1.370 

10 

97.9 

3.2 

48.2 

29.6 

50.6 

29 

57 

50 

85 

135 

1.894 

-1.091 

11 

140.5 

5.2 

52.9 

43.8 

72.9 

28 

46 

41 

100 

141 

3.300 

-1.383 

12 

74.6 

2.2 

33.8 

23.8 

38.4 

22 

63 

60 

74 

134 

1.207 

-0.941 

13 

87.9 

5.0 

38.5 

28.1 

46.5 

20 

54 

47 

96 

143 

1.764 

-0.888 

14 

64.2 

2.7 

32.1 

22.3 

33.5 

23 

72 

49 

89 

138 

1.255 

-0.704 

15 

105.4 

1.5 

44.0 

30.5 

53.5 

27 

50 

49 

71 

120 

2.120 

-1.406 

16 

78.1 

5.6 

40.3 

25.6 

41.9 

24 

59 

57 

65 

122 

1.272 

-1.148 

17 

119.2 

3.5 

57.2 

36.9 

61.4 

29 

57 

43 

82 

125 

2.691 

-1.360 

18 

151.8 

7.7 

74.0 

50.5 

79.8 

27 

51 

39 

83 

122 

3.695 

-1.788 

19 

201.3 

3.4 

90.5 

67.6 

102.4 

23 

54 

47 

79 

126 

4.211 

—2.427 

20 

110.6 

9.6 

59.9 

34.8 

60.1 

23 

74 

49 

91 

140 

2.061 

-1.081 

mean 8 _ 20 

116.2 

5.2 

53.4 

35.9 

60.7 

26.3 

56.1 

48.2 

83.5 

131.6 

2.398 

-1.337 

s 8— 20 

36.7 

2.7 

16.2 

- 

18.8 

5.8 

9.2 

6.2 

10.0 

10.0 

0.958 

0.441 


value equal to approximately 116.2 ± 36.7. Statistically, _correct prediction would be expected 
about 68 percent of the time. A 2-sigma confidence, or R^-^ = 116.2 ± 73.4, should result in 

being correct about 95 percent of the time. Therefore, it would be extremely useful to deduce 
a relationship which gives better estimates for with reduced 1-sigma spreads. In 

APPROACH 1, Rjyj^x values are "detrended" for cycles 8 through 20 by means of standard 

linear-regression time series analysis (e.g., Downie and Heath [168] and Longley-Cook [169]). 

Examination of R^^ values for cycles 8 through 20 clearly shows a downward trend 

between cycles 8 through about 14 and an upward trend from cycle 14 through at least cycle 
19. The entire period spanning cycles 8 through 20 shows an upward trend. Based purely on the 
deduced regression equation, can be "better" predicted since s has been reduced to 33.2. 

Various modifications will be made on this linear-regression fit (denoted LR) to further reduce 
s. These modifications include assuming the presence of an 8-cycle periodicity or an 11-cycle 
periodicity, an 11-cycle periodicity with an imbedded 8-cycle periodicity, and an 11-cycle 
periodicity with an imbedded 2-cycle periodicity. The 8-cycle periodicity was selected because 
it is equivalent to the so-called "Gleissberg" 90-year period often said to be imbedded in the 
sunspot data [55] . The 11-cycle periodicity is the "best-fit" periodicity for the given R^^ 

values (cycles 8 through 20). (The term "best-fit" is used extensively in this report. By this we 
mean that the fit to the observed values of R^^- has the least sum of the squares of the resid- 
uals; see Section in.) The 2-cycle periodicity was selected since throughout much of the period 
spanning cycles 8 through 20, there is an apparent "up-down-up" signal in the data. 
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In APPROACH $, v§iu§§ a?§ fitted for ey§l§s Q through IQ using a gine 

curve with a period equal to 50 cycles. This period was chosen based on the assump- 
tion that the Maunder Minimum actually occurred (the middle of the minimum corres- 
ponds to SCN -6) and that the sine curve is at phase 270 deg at that point and at 
phase 90 deg at SCN 19; hence, the half-period of the sine curve is 25-cycles, 
yielding the period to be 50 cycles. Such an analysis gives a 1-sigma value of 37.0 
and yields a line that is very close in value to the linear-regression line. Modifica- 
tions to this approach are as above in APPROACH 1, using various combinations of 
an 11-cycle, 8-cycle, and/or 2-cycle periodicity imbedded in the 50-cycle fit. 

In APPROACH 3, the LR fit is approximated by a sine curve of period equal to 
90 cycles, such that the curve is at phase 90 deg during the Grand Maximum (at 
SCN -53) and at phase 270 deg during the Maunder Minimum (at SCN -8). Such an 
analysis gives a 1-sigma value of 36.1 and yields a line very close to that of the 
regression line. Modifications to this approach follow that above for APPROACH 1 
and APPROACH 2, using various combinations of an 11-cycle, 8-cycle, and/or 2-cycle 
periodicity imbedded in the 90-cycle fit. 


m- mms mb Bispp§siPN 


Following thg aR£r S gs£ps ffi t0 th ®, 

observed R^ AX ya^jgg fgp pyclgs 8 through 2Q a*| § as *8? fjpst-order 

wed j&gr&g i arke ^d n t ™| 1 i|* m\ je&ft f 1 *’ \mm- 

immm W-Pi8 le m m fit: Isrsisri m irai feetew & n S w 
wtm 1 m te Im 

Mean Fit 


R MAX 116,2 


(1) 


Linear-Regression Fit 

R MA x = 88.1 + 2.0 X .(2) 

50-Cycle Fit 

R max = 100 * 7 + 18 * 7 sin e l (X) (3) 

90-Cycle Fit 

R max = 116 -° + 24 * 7 sin e 5 (X) (4) 


In equation (2), X is the sunspot cycle number (SCN). In equations (3) and (4), 
6^( X) is the angle in degrees for a particular X assuming a 50-cycle periodicity, and 
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0_(X) is the angle in degrees for a particular X assuming a 90-cycle periodicity, 
o 

respectively. These angular values are identified in Table 2. The standard devia- 
tions for each of these fits are 36.7, 36.1, 37.0, and 36.1, respectively. Equation 
(3) is of the form y = a Q + a^ sine, where ap = 100.7 and a^ = 18.7. The constant 

a Q was chosen to be 100.7 on the basis that SCN 19 has the highest known Rjyj AX 
value (= 201.3) and that the lowest possible value for any cycle is 0; thus, a^ is the 
rounded-off half-way point between these two extremes. The constant a^ is the 


"best'-fit" value based on the given a^ and the values of e. (X) as shown in Table 2 

for cycles 8 through 20. The ’’best-fit" a^ was deduced by arbitrarily selecting three 

values of a^ and computing R^^ values for cycles 8 through 20 with the given ap, 

summing the squares of the residuals (where a residual is the difference between the 
observed and computed value of Rjyj AX for a cycle) , and fitting these sums of squares 

of the residuals versus a^ which yields the least value for the sums of the squares 
of the residuals. 


TABLE 2. ANGULAR VALUES (PHASES) FOR SELECTED PERIODIC 
CYCLES AS A FUNCTION OF SCN 


SCN 

CYCLE MAXIMUM 
APPROXIMATE 
YEAR DATE 

o 

> 

o 

1 

s 

^X 

O 

> 

o 

1 

T- 

«r- 

e 2 x 

o 

> 

‘f 

00 

9 3 X 

o 

> 

'i 

CM 

0 4 X 

o 

>• 

0 

1 

o 

o> 

NOTES 

-60 

1088 

241.2 

24.5 

135 

270 

62 



-59 

1099 

248.4 

57.3 

180 

90 

66 



-58 

1110 

255.6 

90.0 

225 

270 

70 



-57 

1121 

262.8 

122.7 

270 

90 

74 



p (APPROXIMATE) 

—56 

1132 

270.0 

155.5 

315 

270 

78 



-55 

1143 

277.2 

188.2 

360 

90 

82 



-54 

1154 

284.4 

220.9 

45 

270 

86 



-53 

1165 

291.6 

253.6 

90 

90 

90 



-52 

1176 

298.8 

286.4 

135 

270 

94 


GRAND MAXIMUM 

—51 

1187 

306.0 

319.1 

180 

90 

98 



-50 

1198 

313.2 

351.8 

225 

270 

102 



-49 

1209 

320.4 

24.5 

270 

90 

106 



-48 

1220 

327.6 

57.3 

315 

270 

no 



-47 

1231 

334.8 

90.0 

360 

90 

114 



-46 

1242 

342.0 

122.7 

45 

270 

118 



-45 

1253 

349.2 

155.5 

90 

90 

122 



-44 

1264 

356.4 

188.2 

135 

270 

126 



-43 

1275 

3.6 

220.9 

180 

90 

130 



-42 

1286 

10.8 

253.6 

225 

270 

134 


(APPROXIMATE) 

-41 

1297 

18.0 

286.4 

270 

90 

138 



-40 

1308 

25.2 

319.1 

315 

270 

142 



-39 

1319 

32.4 

351.8 

360 

90 

146 



-38 

1330 

39.6 

24.5 

45 

270 

150 



-37 

1341 

46.8 

57.3 

90 

90 

154 



-36 

1352 

54.0 

90.0 

135 

270 

158 



-35 

1363 

61.2 

122.7 

180 

90 

162 



-34 

1374 

68.4 

155.5 

225 

270 

166 



-33 

1385 

75.6 

188.2 

270 

90 

170 



-32 

1396 

82.8 

220.9 

315 

270 

174 



-31 

1407 

90.0 

253.6 

360 

90 

178 



-30 

1418 

97.2 

286.4 

45 

270 

182 


(APPROXIMATE) 

-29 

1429 

104.4 

319.1 

90 

90 

186 



-28 

1440 

111.6 

351.8 

135 

270 

190 



-27 

1451 

1 18.8 

24.5 

180 

90 

194 



-26 

1462 

126.0 

57.3 

225 

270 

198 


SPORER 

-25 

1473 

133.2 

90.0 

270 

90 

202 


MINIMUM 

-24 

1484 

140.4 

127.7 

315 

270 

206 



-23 

1495 

147.6 

155.5 

360 

90 

210 



-22 

1506 

154.8 

188.2 

45 

270 

214 



-21 

1517 

162.0 

220.9 

90 

90 

218 


(APPROXIMATE) 

-20 

1528 

169.2 

253.6 

135 

270 

222 


-19 

1539 

176.4 

286.4 

180 

90 

226 



-18 

1550 

183.6 

319.1 

225 

270 

230 



-17 

1561 

190.8 

351.8 

270 

90 

234 



-16 

1572 

198.0 

24.5 

315 

270 

238 
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TABLE 2. (Concluded) 




UJ 

UJ 

UJ 

UJ 

Ul 





-j 

_j 




o 

a 

o 

o 

o 



> 

> 

> 

> 

>- 




a 

o 

a 

o 


CYCLE MAXIMUM 

i 

o 

i 

T— 

1 

00 

CM 

1 

o 


APPROXIMATE 

LO 

*“ 



o 

SCN 

YEAR DATE 

i^X 

0 2 X 

0 3 X 

0 4 X 

05 X NOTES 

-15 

1583 

205.2 

57.3 

360 

90 

242 

-14 

1594 

212.4 

90.0 

45 

270 

246 

-13 

1605 

219.6 

122.7 

90 

90 

250 

-12 

1616 

226.8 

155.5 

135 

270 

254 

-11 

1627 

234.0 

188.2 

180 

90 

258 

-10 

1638 

241.2 

220.9 

225 

270 

262 

_ (APPROXIMATE) 

- 9 

1649 

248.4 

253.6 

270 

90 

266 


- 8 

1660 

255.6 

286.4 

315 

270 

270 


- 7 

1671 

262.8 

319.1 

360 

90 

274 

MAUNDER 

- 6 

1682 

270.0 

351.8 

45 

270 

278 

MINIMUM 

- 5 

1693 

277.2 

24.5 

90 

90 

282 


- 4 

1704 

284.4 

57.3 

135 

270 

286 


- 3 

1715 

291.6 

90.0 

180 

90 

290 J 

_ (APPROXIMATE) 

- 2 

1726 

298.8 

122.7 

225 

270 

294 

- 1 

1737 

306.0 

155.5 

270 

90 

298 

0 

1748 

313.2 

188.2 

315 

270 

302 

1 

1759 

320.4 

220.9 

360 

90 

306 

2 

1770 

327.6 

253.6 

45 

270 

310 

3 

1781 

334.8 

286.4 

90 

90 

314 

4 

1792 

342.0 

319.1 

135 

270 

318 

5 

1803 

349.2 

351.8 

180 

90 

322 

6 

1814 

356.4 

24.5 

225 

270 

326 

7 

1825 

3.6 

57.3 

270 

90 

330 

8 

1836 

10.8 

90.0 

315 

270 

334 

r 

9 

1847 

18.0 

122.7 

360 

90 

338 

1 BEST-QUALITY 

10 

1858 

25.2 

155.5 

45 

270 

342 

T SUNSPOT DATA 

11 

1869 

32.4 

188.2 

90 

90 

346 

12 

1880 

39.6 

220.9 

135 

270 

350 

13 

1891 

46.8 

253.6 

180 

90 

354 

14 

1902 

54.0 

286.4 

225 

270 

358 

15 

1913 

61.2 

319.1 

270 

90 

2 

16 

1924 

68.4 

351.8 

315 

270 

6 

17 

1935 

75.6 

24.5 

360 

90 

10 

18 

1946 

82.8 

57.3 

45 

270 

14 

19 

1957 

90.0 

90.0 

90 

90 

18 

20 

1968 

97.2 

122.7 

135 

270 

22 

21 

1979 

104.4 

155.5 

180 

90 

26 

22 

1990 

111.6 

188.2 

225 

270 

30 

23 

2001 

118.8 

220.9 

270 

90 

34 

24 

2012 

126.0 

253.6 

315 

270 

38 

25 

2023 

133.2 

286.4 

360 

90 

42 

26 

2034 

140.4 

319.1 

45 

270 

46 

27 

2045 

147.6 

351.8 

90 

90 

50 

28 

2056 

154.8 

24.5 

135 

270 

54 

29 

2067 

162.0 

57.3 

180 

90 

58 

30 

2078 

169.2 

90.0 

225 

270 

62 


Similarly, equation (4) is of the form y = a Q + a.^ sine, where a Q = 116 and 

= 24.7. The constant ag was selected to be approximately the value that the 

linear-regression equation [eq. (2)] yields for SCN 14, the mid-point of the data set. 
As above, a^ is the "best-fit" for the given Sq and the values of 6 (X) shown in 

Table 2, using the least-squares method. 

In Figure 2, the residuals are plotted in units of standard deviation for each of 
the fits shown in Figure 1. For convenience, the figure is arbitrarily divided into 
several smaller regions. The middle region of each plot is designated "average" 
since approximately 50 percent, on average, of the R MAX values are expected to lie 

within its bounds; the region just above (and below) "average" is called "above (below) 
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Figure 2. Statistical spread versus SCN (cycles 8 through 20), 














average and approximately 20 percent, on average, of the cycles are expected to fall 
within its bounds. The top-most (bottom-most) region is called "anomalous high (low) 
and represents the approximately upper (lower) 5 percent, on average, of the R 

1V1AA 


values which are greater than 1.645 standard deviations above (below) the computed 
value for R MAX - Based on these single-component fits, SCN 19 is observed to 


Tf 


probably be an anomalously high R MAX 


cycle * 


In subsections A, B, and C which follow, the effects of selected modifications 
to the aforementioned basic curve-fits (linear regression, 50-cycle, and 90-cycle fits) 
are examined. In subsection D, computed values of Rjyj AX are compared using the 


preferred empirical fit with values reported for earlier epochs (SCN 1 through 7 and 
the Maunder Minimum). (The reader is reminded that phase for the empirical fit is 
based, in part, on the occurrence of the Maunder Minimum; however, amplitude for 
the fit is based only on_cycles 8 through 20.) In Subsection E, the empirical fit is 
used to infer values of Rjyj AX for future epochs (SCN 21, 22, and beyond). 


A. Linear-Regression Fits (APPROACH 1) 

Following the first modification in APPROACH 1, Figure 3 shows the linear- 
regression fit with an 8-cycle variation. The equation for this fit is: 


R MAX = 88,1 + 2 -° x + 19 - 8 sin V x > 


(5) 
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Figure 3. LR + 8 FIT. 
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where X is again SCN and 0g(X) is found in Table 2. The constant 19.8 is the 


"best-fit" for the values given in the table and above. This fit only slightly reduces 
the statistical spread (s = 33.4 for the LR + 8 FIT as compared to ej = 36.1 for the 
LR FIT alone); so, if indeed there is an 8-cycle periodicity in the Rj^x recor< l for 

cycles 8 through 20, it does not appear to contribute much to the variation in R MAX 


value with time. Based on this two-component fit, both cycles 8 and 19 are observed 
to probably be anomalously high. 


A better two-component fit is plotted in Figure 4 which shows the linear- 
regression fit with an 11-cycle variation. Its equation is: 


R MAX = 88,1 + 2 -° x + 30.7 sin e 2 (x > 


( 6 ) 


As before, X is SCN and 0g(X) can be found in Table 2. The constant 30.7 is the 

"best-fit" for the given table values and constants used in equation (6). This fit 
reduces the statistical spread by about 30 percent (s = 25.2 for the LR + 11 FIT). 
Based on this fit, cycle 19 is anomalously high and cycle 20 is anomalously low. 

Figure 5 shows a three-component fit based on the linear-regression fit with an 
11-cycle and 8-cycle variation. Its equation is: 


R MAX = S8.1 + 2.0 X + 30.7 sin 0 2 (X) + 11.4 sin 0g(X) 


(7) 
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Figure 4. LR + 11 FIT. 
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Figure 5. LR + 11 + 8 FIT. 

Again, this equation only slightly reduces the statistical spread as compared to equa- 
tion (6) (s = 24.9 for the LR + 11 + 8 FIT). Based on this fit, cycle 20 is anomal- 
ously low. 

The final fit following APPROACH 1 is shown in Figure 6. It is a three- 
component fit based on the linear-regression fit with an 11-cycle and 2-cycle varia- 
tion. Its equation is: 


R M ax = 88.1 + 2.0 X + 30.7 sin 0 2 (X) + 15.8 sin 8 4 (X) 


( 8 ) 


This fit reduces the original statistical spread from s = 36.1 to s = 19.7, a 45-percent 
reduction. (The variance has been reduced by 70 percent.) No cycles are observed 
to be anomalously low or high based on this fit, although several are in the above 
average or below average categories. 


B. 50-Cycle Fits (APPROACH 2) 

In APPROACH 2, the baseline fit is a 50-cycle periodicity. Modifications to it 
will be made as in APPROACH 1. Figure 7 depicts an 8-cycle variation on the 50- 
cycle fit. Its equation is: 


R MAX = 100,7 + 18,7 sin 9 1 (X) + 21 sin e 3 (X) 


(9) 
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Figure 7. 50+8 FIT. 










The statistical spread about this line is only slightly better than the 50-cycle fit 
alone (s = 33.8 for the 50 + 8 FIT). Cycles 8 and 19 both appear as anomalously 
high cycles . 

In Figure 8, an 11-cycle variation on the 50-cycle fit is plotted. As in 
APPROACH 1, a much better fit is suggested. Its equation is: 


R MAX = 100,7 + 18 * 7 sln 8 i( x ) + 33.2 sin 0 2 (X) . (10) 


Its statistical spread is s = 24.5, about 34 percent more restrictive than the 50-cycle 
fit alone. Cycle 19 is the only cycle that appears to be anomalously high. 
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Figure 8. 50 + 11 FIT. 

In Figure 9, an 11-cycle and an 8-cycle variation have been placed on the 
50-cycle fit. Its equation is: 


R MAX = 100,7 + 18,7 sin 0 1 (X) + 33,2 sin e 2 (X) + 8,7 sin e 3< x > . (11) 


As before, the presence of an 8-cycle periodicity has little effect upon the reduction 
of the statistical spread (s = 23.8 for the 50+11+8 FIT). Cycles 19 and 20 appear 
to be anomalously high and low, respectively. 
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Figure 9. 50+11+8 FIT. 

In Figure 10, an 8-cycle and 2-cycle variation has been included with the 
50-cycle periodicity. Its equation is: 


R MAX = 100 - 7 + 18 - 7 sin 6j(X) + 21*9 sin 0 3 (X) + 12.8 sin e 4 (X) . (12) 

This fit is worse than the preceding fit, having s = 31.2. Cycle 8 appears to be an 
anomalously high cycle. 

A much better fit is that shown in Figure 11 where an 11-cyele and a 2-cycle 
variation has been added to the baseline 50-cycle fit. Its equation is: 


R MAX = 100 * 7 + 18 ’ 7 sin 0 1 (X) + 33,2 sin 0 2 (X) + 15,4 sin 0 4 (X) * (13) 


Its statistical spread is 49 percent smaller than that deduced for the baseline (s = 
19.0 for the 50+11+2 FIT; the variance has been reduced by 74 percent). Only 
cycle 19 appears to be anomalously high. 


C. 90-Cycle Fits (APPROACH 3) 

In this subsection modifications to a 90-cycle periodicity, which approximates 
the linear-regression equation used in APPROACH 1, are examined. Figure 12 shows 
an 8-cycle variation on the 90-cycle fit. Its equation is: 
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Figure 12. 90 + 8 FIT. 

R MAX = 113-0 + 24.7 s ^ n e 5^ X ^ + 23 *3 s ^ n 0 3^ x ^ • (14) 

Its spread is 33.1, only slightly reduced from the 90-cycle baseline spread of 36.1. 
Cycles 8 and 19 appear to be anomalously high. 

A better fit is obtained using an 11-cycle variation on the 90-cycle periodicity, 
shown in Figure 13. Its equation is: 



r max * H6.0 + 24.7 sin 0g(X) + 35.4 sin 0 2 (X) . (15) 

Its 1-sigma spread is 24.6 units. Cycles 19 and 20 appear to be anomalously high 
and low, respectively. 

Figure 14 shows an 8-cycle and 2-cycle variation on the 90-cycle periodicity. 
Its equation is 


R MAX = 116,0 + 24,7 sin e 5 (X) + 20,3 sin e 3 (X) + 13,2 sin 0 4 (X) ' (16) 

Its 1-sigma spread is 30.4 and cycle 8 appears to be anomalously high. 
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Figure 13. 90 + 11 FIT. 
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Figure 14. 90+8+2 FIT . 
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Figure 15 depicts the function which had the smallest spread for the examples 
tested. It is an 11-cycle and 2-cycle variation on the 90-cycle periodicity. Its 
equation is: 


R MAX = 116 *° + 24,7 e 5 ^ x ^ + 35.4 sin 6 2 (X) + 15.9 sin 0 4 (X) * (17) 


Its 1-sigma spread is 18.8 units, a 48-percent reduction from the baseline spread. 
(The 'variance has been reduced by 73 percent.) Cycle 17 appears to be anomalously 
low. 
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Figure 15. 90 + 11 + 2 FIT. 


D. Pre-Cycle 8 Fits 

According to the above analysis, values of R M ^x ^ or c Y cles 8 through 20, the 

most reliably known, can be fit very well by equation (17). Using this preferred 
empirical function, one can now extrapolate ’'backwards” in time to determine computed 
smoothed sunspot number at cycle maximum for earlier epochs, (It is realized that 
equation (17) is merely a fit for a small sample, _and that this fit may have no ”post- 
dictive." or "pre-dictive" value.) In Figure 16, Rjyj^x * s com puted backwards in time 

to about A.D. 1080, corresponding to SCN -60. The top-most line is the 90+11+2 
FIT and the bottom-most line is the 90 + 8 + 2 FIT, drawn for comparison. Observed 
values of R MAX for cycle 1 through 21 are drawn as open circles; computed values 
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are drawn as darkened triangles. The approximate dates, as stated by Eddy [3,4], 
for the Maunder Minimum, Sporer Minimum, and Grand Maximum are shown. A table 
legend is included which identifies some major milestones in solar research, as related 
to the study of sunspots. Computed mean values and standard deviations are shown 
for the two minima and single maximum periods. A note reminding us that the best- 
quality sunspot data extends from about cycle 8 to the present is also included. One 
quickly observes that, while the 90+11+2 FIT gives a very good representation of 
cycles 8 through 21, it gives a very poor match to cycles 2 through 7, sometimes 
departing by more than 3 sigma from computed values. Cycles 2 through 7 appear 
to be better- fitted by the 90 + 8+2 FIT . This seems to be supportive of the belief 
that there is an 8-cycle periodicity imbedded in the sunspot record; however, how 
much weight should be given to this view and even to the firmness of the R MAX 

values for cycles 2 through 7 remain unclear . Additional investigation should be 
pursued . 

Table 3 gives values for several selected cycle-related parameters for cycles 1 
through 7. Mean values and standard deviations are given for the group cycles 1 
through 7 and for the groups cycles 8 through 20 and 1 through 20. Clearly, when 
parametric values are compared for the different groupings, it is found that in some 
cases they are quite dissimilar; e.g., ASC and DES. This may be an indication that 
smoothed sunspot number at cycle maximum for cycles observed prior to cycle 8 are 
in error. (It is known from Waldmeier [170] that daily records of sunspot number, 
which form the basis for a determination of monthly mean sunspot number which is 
then used to compute smoothed sunspot number, are incomplete.) Rather, it may 
mean that the period encompassing cycles 2 through 7 was an anomalous timespan, 
with values departing from computed values by substantial amounts. (Very recently, 
Sonett [36] has addressed this particularly anomalous era and found the sunspot 
record to be consistent with the idea that, perhaps, there may have been a "book- 
keeping" error in the compilation of sunspot number about 1780 to 1800.) 


TABLE 3. SELECTED PARAMETRIC VALUES FOR CYCLES 1 THROUGH 7 


SC N 

MINIMUM 

OCCURRENCE DATE 

MAXIMUM 
OCCURRENCE DATE 

r min 

r max 

ASC 

DES 

MIN— MIN 
PERIOD 

1 

FEB 1755 

JUN 1761 

8.4 

86.5 

75 

60 

135 

2 

MAY 1766 

SEP 1769 

11.2 

115.8 

39 

69 

108 

3 

JUN 1775 

MAY 1778 

7.2 

158.5 

35 

76 

111 

4 

SEP 1784 

FEB 1788 

9.5 

141.2 

41 

96 

137 

5 

APR 1798 

FEB 1805 

3.2 

49.2 

81 

65 

146 

6 

JUL 1810 

MAY 1816 

0.0 

48.7 

70 

83 

153 

7 

APR 1823 

NOV 1829 

0.1 

71.7 

79 

48 

127 


MEAN“j_7 

5.7 

95.9 

60.0 

71.0 

131.0 

s 1— 7 

4.2 

40.4 

19.1 

14.6 

15.6 

MEAN 8 _2Q 

5.2 

116.2 

48.2 

83.5 

131.6 

o 

CNI 

1 

CO 

</> 

2.7 

36.7 

6.2 

10.0 

10.0 

MEAN*j„20 

5.4 

109.1 

52.3 

79.1 

131.4 

s 1 —20 

3.3 

39.2 

13.6 

13.2 

12.2 
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Of the three aforementioned, preceding epochs, only the Maunder Minimum can 
be examined more closely. According to equation (17), as Figure 16 shows, a minimum 
R MAX value occurs i n the 17th century. However, the values postdicted by equation 

(17) are higher than has been suggested for this period (e.g., Eddy [1,8]). (Equa- 
tion (17) suggests a mean value for R MAX of about 91 ± 32.) Recently, Xu and 

Jiang [35] re-assessed the 17th Century Chinese sunspot records and concluded that 
the inference to a Maunder Minimum may be incorrect (Link [23] and Schove [24]). 
They found sightings which allowed them to argue that solar cycles were continuing 
during the Maunder Minimum and that the average annual sunspot number at cycle 
maximum during this period was around 50 or less. Very recently, Zhang [37], 
based on Korean auroral records for the period A.D. 1507-1747, reports that solar 
activity during the Maunder Minimum may not have been as low as Eddy has claimed. 
Also, Hameed and Wyant [38] have found an approximate 23-year cycle in surface 
temperatures in central England during and after the Maunder Minimum, suggesting 
that if variation in climate may be partly due to solar activity, then solar activity 
may still have been very much apparent even during the Maunder Minimum. Figure 
16 gives R max values which can be converted to annual numbers for direct compari- 
son. Figure 17 plots annual sunspot number at cycle maximum, R(A), versus R]y[ AX > 



r max, smoothed sunspot number at cycle maximum 


Figure 17. R(A) versus R]y[ AX - 
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based on cycles 8 through 20. The Pearson correlation coefficient r is 1.00 with a 
standard error of estimate S equal to about 3.8. The equation relating the two 
parameters is : y 


R(A) = -1.063 + 0.952 R MAX 


( 18 ) 


Based on equation (18) and the 90 + 11 + 2 FIT, R(A) values are expected to have 
been -about 86 + 31 during the Maunder Minimum. If the Maunder Minimum is viewed 
as an anomalous timespan where R(A) values are depressed about 1 or more standard 
deviations below normal, then R(A) is expected to often be below 55, very near 
values reported by Xu and Jiang. 


In Table 4, the Xu and Jiang results are recast along with those reported by 
Schove [22] and R(A) estimates based on the 90+11+2 and 50 + 11 + 2 FITS. 

The table spans SCN -14 to SCN -3, where the Maunder Minimum approximately spans 
SCN -9 to SCN -3. Six of the ten records (Xu and Jiang only give results for SCN 
-14 to SCN -5) for which Xu and Jiang give results are very compatible with those 
obtained using either the 90 + 11+2 or 50 + 11 + 2 FITS. If the statistical spread 
is allowed to be greater than 1-sigma, then, essentially all the records overlap. 

This comparison suggests that while a Maunder Minimum very_ well may have occurred, 
solar periodicities were probably still working and values of R MAX given by the 


empirical function (depressed by at least 1-sigma, possibly more) may be representa- 
tive of the period. (Recently, Eddy [8] has re-appraised the Maunder Minimum and 
argues that the basis for the occurrence of such a minimum remains firm, with peak 
annual sunspot number being of the order of 20 or below during this protracted 
period of low sunspot number.) 


TABLE 4. COMPARISON OF REPORTED AND COMPUTED R . VALUES 

FOR THE MAUNDER MINIMUM MAA 


SCN 


R' A » 2 est 

r(a)3 comp 

r(a)4 comp 

NOTES 

-14 

50 

80 

106 

102 



-1:3 

>160 

90 

129 

125 



-12 

85 

100 

85 

80 



-11 

160 

70 

97 

91 



-10 

50 

40 

48 

44 



- 9 

50 

j.50 

68 

63 



- 8 

45 

1 

38 

33 



- 7 

<45 

60 

79 

71 



- 6 

45 

50 

66 

58 

MAUNDER 






MINK 

\/IUM 

- 5 

<45 

30 

116 

105 



- 4 

* 

50 

100 

89 



- 3 

* 

130 

135 

125 

’ 

L_ 


*No value given in Xu and Jiang [35] 

1 According to Xu and Jiang [35], based on naked-eye Chinese 
sunspot records 
^According to Schove [22] 

^Based on 90 + 11 + 2 FIT 
4 Based on 50 + 11 + 2 FIT 
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E. Cycles Beyond SCN 20 


This report will now extrapolate the empirical function "forwards" in time, at 
least for cycles 21 and 22, the two cycles of more immediate concern. The 90 + 11 + 

2 FIT predicts that cycle 21 has a 1-sigma confidence R MAX value equal to 157.3 ± 

18.8. This value is remarkably consistent with that predicted by Sargent [64] (= 154) 
and is very close to that which has now been observed (= 164.5), Table 5 lists 


estimates for R 


MAX 


for cycle 21 based on all the approaches given earlier; it also 


gives residuals. It was seen that cycle 21 was about 0.38 s unit above the predicted 
value, in the average range of Figure 15. Table 6 lists the cycle-related parameters 
projected for cycle 21, based on the 90 +11 + 2 FIT and on the actually observed 
R max value. (Values for the cycle-related parameters are based on regression equa- 


tions using R max 


as the independent variable, given in Figures 18 through 20, 


included_in the Appendix.) Also given in Table 6 are the occurrence dates for 
Rmax’ ^CHM ( asceirt an( l descent side), and cycle end date, denoted here as 

R M i N (SCN + 1) occurrence which marks the beginning of the next cycle (SCN 22). 


TABLE 5. R max (COMP) VALUES FOR SCN 21 


METHOD R MAX (COMP) ± 1 SIGMA RES = R MAX (OBS) - R MAX (COMP) 


MEAN 

116.2 ±36.7 

48.3 

LR 

130.1 ± 36.1 

34.4 

LR + 8 

130.1 ± 33.4 

34.4 

LR + 11 

142.8 ± 25.2 

21.7 

LR + 11+8 

142.8 ± 24.9 

21.7 

LR + 11 + 2 

158.6 ± 19.7 

5.9 

50 

118.8 + 37.0 

45.7 

50 + 8 

118.8 ± 33.8 

45.7 

50 + 11 

132.6 ± 24.5 

31.9 

50+11+8 

132.6 ± 23.8 

31.9 

50 +8 + 2 

131.6 ±31.2 

32.9 

50 + 11+2 

148.0 ± 19.0 

16.5 

90 

126.8 ± 36.1 

37.7 

90+ 8 

126.8 + 33.1 

37.7 

90 + 11 

141.5 ± 24.6 

23.0 

90+ 8 + 2 

140.0 ± 30.4 

24.5 

90 + 11+2 

157.3 ± 18.8 

7.2 


Table 6 reveals that, based on the 90 +11+2 FIT R jyj AX projection (= 157.3), 

eycle_21 maximum would have been estimated to have occurred about February 1980 
and occurrence dates about October 1978 and January 1983; cycle end date 


Based on actual dates of occurrence, 

it is noted that R n) , AV occurrence would have been missed by only 2 months (December 


would be expected to occur about March 1987. 
l MAX 


1979), and K 


CHM 


ascent-side and descent-side occurrences by 4 months and 2 months, 


respectively (June 1978 and March 1983). Cycle 21 is still in decline, so a com- 
parison between computed and observed cycle end dates is not yet possible. Wilson 
[44] and Wilson et al. [45,100] have discussed the expected cycle duration for 
cycle 21. If cycle 21 is a short-period cycle (cycle duration less than 128 months). 
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TABLE 6. SELECTED PARAMETRIC VALUES FOR SCN 21 


PARAMETER 

MEAN 

VALUE 1 

STANDARD 

DEVIATION 1 

OBS 

VALUE 

COMP. 

VALUE 2 

COMP. 

VALUE 3 

^MAX 

116.2 

36.7 

164.5 

157.3 


S MIN 

5.2 

2.7 

12.2 

6.1 

6.2 

Rmean 

53.4 

16.2 


70.7 

73.7 

S MEAN 

35.9 

12.2 


48.3 

50.4 

S CHM 

60.7 

18.8 

88.4 

81.7 

85.4 

asc chm 

26.3 

5.8 

24 

27.6 

27.8 

d chm 

56.1 

9.2 


50.9 

50.0 

ASC 

48.2 

6.2 

42 

44.3 

43.6 

DES 

83.5 

10.0 


84.5 

84.7 

MIN— MIN 






PERIOD 

131.6 

10.0 


128.8 

128.3 

SLOPEasc 

2.398 

0.958 

3.626 

3.426 

3.606 

slope des 

-1.337 

0.441 


-1.767 

-1.851 

R MAX OCCURRENCE DATE: 

R CH M OCCURRENCE DATES: 

R M|N (SCN + 1) OCCURRENCE DATE: 

ASC-SIDE 

DES-SIDE 

QB§, 

DEC 1979 

JUN 1978 
MAR 1983 

COMP. 2 
FEB 1980 

OCT 1978 
JAN 1983 
MAR 1987 

COMP. 3 

OCT 1978 
DEC 1982 
FEB 1987 


1 BASED ON CYCLES 8 THROUGH 20 

2 BASED ON 90+11+2 FIT; R max (COMP> = 157.3 

3 BASED ON R max (OBS) = 164.5 


then cycle 21 should end before March 1987 (best-guess date: August 1986); other- 

wise, cycle 21 will be a long-period cycle (cycle duration greater than 133 months) 
ending no earlier than August 1987 (best-guess date: February 1988). 

In Table 7, computed R MAX values with 1-sigma confidence for cycle 22 are 

listed. Extrapolation of the 90 + 11+2 FIT yields an estimate of R MAX for cycle 22 

to be 107.4. (Sargent and Heckman [66] and Sargent [171] have discussed current 
efforts to predict cycle 22; also, a Solar Terrestrial Prediction Workshop with the 
expressed objective of discussing current efforts to predict cycle 22 occurred 18-22 
June 1984 in Meudon, France.) Using this value, Table 8 projects values for the 
cycle-related parameters for cycle 22. The R^^ occurrence date is projected to be 
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TABLE 7. R max (COMP) VALUES FOR SCN 22 


METHOD 

R MAX (C0MP > 1 1 SIGMA 

LR + 1 1 + 2 FIT 

111.9 ± 19.7 

50 + 11+2 FIT 

98.0 + 19.0 

90+ 11 +2 FIT 

107.4 ± 18.8 

90+ 8 + 2 F IT 

100.8 ± 24.5 


TABLE 8. 

SELECTED 

PARAMETRIC VALUES FOR 

SCN 22 

PARAMETER MEAN VALUE 1 

STANDARD DEVIATION 1 

COMPUTED VALUE 

^MAX 

116.2 

36.7 

107.4 

^MIN 

5.2 

2.7 

5.0 

^MEAN 

53.4 

16.2 

49.7 

S MEAN 

35.9 

12.2 

33.3 

^CHM 

60.7 

18.8 

56.2 

asc chm 

26.3 

5.8 

26.0 

d chm 

56.1 

9.2 

57.2 

ASC 

48.2 

6.2 

49.0 

DES 

83.5 

10.0 

83.2 

MIN— MIN 
PERIOD 

131.6 

10.0 

132.2 

sl° pe asc 

2.398 

0.958 

2.178 

slope des 

—1.337 

0.441 

-1.245 

R msn OCCURRENCE DATE: 


COMP. 
MAR 1987 


RjViax OCCURRENCE DATE: 


APR 1991 


R CHM OCCURRENCE DATES: ASC-SIDE 

DES-SIDE 

MAY 1989 
FEB 1994 


R M|N (SCN + 1) OCCURRENCE DATE: 

MAR 1998 



’BASED ON CYCLES 8 THROUGH 20 

2 BASED ON 90 + 1 1 + 2 FIT; R MAX (COMP) = 107.4 
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April 1991, occurrence dates May 1989 and February 1994, and the cycle-end 

date March 1998, all dates based on the March 1987 projected end date for cycle 21. 
Other dates, based on bifurcation of sunspot cycles into short -period and long-period 


cycles , have been summarized elsewhere ; (Wilson a] . 

cycle 21 which is the second highest documented R MA - 


[45,100]). In contrast to 
cycle, cycle 22 is projected 


to be a low-valued R MAX cycle, very probably of long cycle duration. Future cycles 


are projected to be decreasing in R MAX value through at least cycle 24 and then to 


be increasing in R MAX value thereafter, at least through cycle 30 (Fig. 16). The 

next "Grand Maximum" like period is projected to be around SCN 40, corresponding 
to ca. A.D. 2180. 


IV. CONCLUSIONS 


This study has shown that values of Rjyj AX for modern sunspot cycles (i.e., 

cycles 8 through 21) can be satisfactorily fitted using a three-component sine curve 
representation, having an 11-cycle and a 2-cycle variation imbedded in an approxi- 
mately 90-cycle periodicity (in terms of sunspot cycle number) . The 90+11 + 2 FIT 
yields "minima," non- zero, low-value R MAX periods, such as the Maunder and Sporer 

Minima, and a "maximum," high-value Rjyj AX period, such as the Grand Maximum, as 

reported by Eddy [3,4], although the values of Rm AX projected for the protracted 

sunspot minimum periods are considerably higher than has been estimated by other 
investigators [8,35]. Annual sunspot numbers, as reported by Xu and Jiang [35], 
for the Maunder Minimum are found to be very close to values "post-dieted" by the 
90 + 11 + 2 FIT, especially if we consider the Maunder Minimum (and other such 
minima) to be a brief period when sunspot number is about 1-sigma below "normal" 
(i.e., below computed values). The 90 + 11 + 2 FIT, extrapolated forward in time, 
has been used to give estimates of Rjyj AX for cycles 21 and 22, and a number of 

useful linear relationships based on Rjyj AX [44] has been applied to project the 

behavior of these cycles. A comparison of observed and predicted values for the 
limited number of known cycle-related parameters for cycle 21 yields very good 
agreement. Cycle 22 is predicted to be a low-valued R... „ /long-period cycle, with 
Rmax ec l ua l to about 107. 
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